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Abstract. We study the asymmetric zero-range process (ZRP) with L sites and 
open boundaries, conditioned to carry an atypical current. Using a generalized Doob 
/i-transform we compute explicitly the transition rates of an effective process for which 
the conditioned dynamics are typical. This effective process is a zero-range process with 
renormalized hopping rates, which are space dependent even when the original rates are 
constant. This leads to non-trivial density profiles in the steady state of the conditioned 
dynamics, and, under generic conditions on the jump rates of the unconditioned ZRP, 
to an intriguing supercritical bulk region where condensates can grow. These results 
provide a microscopic perspective on macroscopic fluctuation theory (MFT) for the 
weakly asymmetric case: It turns out that the predictions of MFT remain valid in the 
non-rigorous limit of finite asymmetry. In addition, the microscopic results yield the 
correct scaling factor for the asymmetry that MFT cannot predict. 
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There has been considerable progress in recent years in the understanding of 
nonequilibrium, current-carrying systems [1-5] . Recently, much effort has been devoted 
to the study of fluctuations which result in atypical currents in such systems [1,5], where 
equilibrium concepts of entropy and free energy could be generalized to a nonequilibrium 
setting. Following the seminal papers [6, 7], the non-typical large-deviation behaviour 
of Markovian stochastic lattice gas models has come under intense scrutiny in the 
framework of what is now known as macroscopic fluctuation theory (MFT) [5]. In 
particular, the asymmetric simple exclusion process (ASEP) [8, 9] has been studied in 
great detail under the condition that the dynamics exhibits a strongly atypical current 
for a very long interval of time. Among the questions of considerable importance is 
the optimal density prohle that realizes such a rare large deviation. Using MFT it was 
found for the weakly asymmetric simple exclusion process, where the hopping bias is 
of order 1/L, that for an atypically low current a dynamical phase transition occurs 
in the case of periodic boundary conditions: Above a critical non-typical current the 
optimal macroscopic density prohle is hat, while rare events below that critical current 
are typically realized by a traveling wave [10]. This phenomenon was subsequently been 
studied numerically with Monte-Carlo simulations [11, 12]. Non-hat optimal prohles 
were obtained also for open boundary conditions [13]. 

More recently, the space-time realizations of large-deviation events in the ASEP 
with finite (strong) hopping bias were studied on the microscopic scale. This was done 
by conditioning the time-integrated current in a finite time interval to attain some 
atypical value. Furthermore, the full conditioned probability distribution on a finite 
lattice has been examined. For low non-typical current, the microscopic structure of 
a huctuating traveling wave and an antishock has been identihed [14-16]. Recently 
duality was used to extend the microscopic approach to show that the traveling wave 
may exhibit a microscopic hne structure consisting of several shocks and antishocks 
[17]. These hndings are consistent with the macroscopic results for the optimal density 
prohle [10, 13], but provide much more information in that the complete time-dependent 
measure was obtained for the microscopic dynamics. 

Using a type of Doob’s h-transform [18-21] one may also study the large time limit 
on microscopic lattice scale. This transformation generates ehective dynamics which 
make the rare large deviation behaviour of the original dynamics typical. For the ASEP 
with periodic boundary conditions and large non-typical current, inaccessible to MFT, it 
was found that a diherent dynamical phase transition involving a change of dynamical 
universality class occurs: Instead of the well-known generic universality class of the 
KPZ equation with dynamical exponent z = 3/2 for the typical dynamics [22] one has 
a ballistic universality class with dynamical exponent z = 1 where huctuations spread 
much faster [23-25]. In fact, it turns out that these results follow predictions from 
conformal field theory and are thus expected to be universal [26]. Moreover, using the 
h-transform technique, long-range effective interactions which make these rare events 
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typical were obtained [23]. 

In this paper, we study similar questions, regarding the emergence of atypical 
currents, in another prototypical model of nonequilibrium systems — the zero range 
process (ZRP) [27, 28]. The steady state of the ZRP is exactly soluble, and thus it has 
often been used in studies of out-of-equilibirum features. The ZRP is a stochastic lattice 
gas model where each lattice site i can be occupied by an arbitrary number > 0 of 
particles. Particles move randomly to neighbouring sites with a rate u{ni) that depends 
only on the occupation number of the departure site and not on the state of the rest 
of the lattice. We consider a one-dimensional lattice of sites i = 1,... ,L with open 
boundaries where the system exchanges particles with external reservoirs. The hopping 
events may be biased in one direction. 

The precise definition of the model is given below. As an introduction we summarize 
some well-known features of the ZRP. The stationary state of the model with periodic 
boundary conditions, where the total particle number is conserved, may exhibit a 
condensation transition above a critical density. In this scenario, the existence of which 
depends on the form of the hopping rates u{n), all sites but one are occupied by a 
particle number that fluctuates around the critical density, while one randomly selected 
site carries all the remaining excess particles, whose number is of the order of the lattice 
size (for a review see [28]). On the other hand, for open boundary conditions where 
the total particle number fluctuates due to random injection and absorption of particles 
at the two boundaries of the system, the scenario is different: For boundary rates for 
which a stationary distribution exists condensation never occurs. Only in the non¬ 
stationary case of very strong injection a dynamical condensation phenomenon occurs: 
The boundary sites can become supercritical, in which case the occupation number of 
these sites grows indehnitely [29]. This gives rise to a non-stationary condensate-like 
structure, but, in further contrast to the usual stationary bulk condensate, this can 
happen for any choice of bounded hopping rates u{n) and, moreover, the phenomenon 
is strictly limited to the boundary sites. 

This brief outline of the properties of the ZRP concerns typical behaviour of the 
stochastic dynamics. It is the purpose of this work to study on microscopic scale 
the behaviour of the ZRP in a regime of strongly atypical behaviour of the particle 
current. This is of interest for several reasons. First, it serves not only as a test for 
macroscopic fluctuation theory (MFT) [5], but also probes its validity in the regime 
of strong asymmetry where MFT is not rigorous. Second, we will be able to compute 
effective microscopic interactions that make atypical behaviour typical. Third, it will 
transpire that unexpected and novel condensation patterns can occur even when the 
unconditioned dynamics do not exhibit condensation. 

Current large deviations of finite duration have been investigated for the ZRP 
in the context of the breakdown of the Gallavotti-Cohen symmetry for the current 
distribution in a ZRP with open boundary conditions [30, 31]. It turns out that the 
failure of the Gallavotti-Cohen symmetry argument, which is based on a very general 
time-reversal property of stochastic dynamics [32, 33], can be related to the formation 
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of “instantaneous condensates” [30]. These condensates were investigated recently in 
terms of Doob’s h-transform for the ZRP with a single site [34]. It was shown that in 
some parameter regimes, Doob’s h-transform fails to represent the effective dynamics 
that make the large deviation typical. 

The microscopic large-deviation properties of the ZRP in the regime of atypical 
currents that persist for a very long time have not yet been explored. In this paper 
we address this problem. Moreover, we do not limit ourselves a single site; rather, 
we consider the full problem of open boundaries with any number of sites. It turns 
out that Doob’s h-transform can be computed exactly from a product ansatz for the 
lowest eigenvector. Thus we are able to calculate the effective interactions that make 
the current large deviations typical. This turns out to be a ZRP with space-dependent 
hopping bias. Interestingly, the effective process satisfies detailed balance if and only of 
one conditions on vanishing macroscopic current. These results are outside the scope of 
MFT. 

Somewhat surprisingly the exact results show that bulk condensation in the 
conditioned open system may occur, as is the case for periodic boundary conditions 
under typical dynamics. However, in contrast to typical behaviour (and to naive 
expectation), the results suggest that a whole lattice segment may become supercritical 
rather than just a single site. The segment location and length are fixed by the current 
on which one conditions rather than randomly fluctuating as one has for the condensate 
position in the periodic case with typical dynamics. 

Our exact results are valid for any asymmetry, including the weakly asymmetric 
case. This allows for a comparison with the predictions of the macroscopic fluctuation 
theory which we apply to the ZRP with weak asymmetry. The microscopic results show 
that the MFT results remain intact in the limit of finite asymmetry, a limit in which 
the validity of the macroscopic approach is a-priori questionable. Moreover, we obtain 
the scaling factor for the asymmetry that cannot be computed from the MFT. 

The paper is organized as follows: In Sec. 2 we introduce the model and define 
the conditioned dynamics. In Sec. 3 we derive the exact microscopic results for the 
h-transform. This yields the effective dynamics and the spatial condensation patterns 
of the conditioned process. In Sec. 4 we follow the macroscopic approach for the weakly 
asymmetric case to compute the optimal macroscopic profile that realizes a current large 
deviation. 

2. Open ZRP conditioned on an atypical cnrrent 
2.1. Definition of the model 

In the bulk of the lattice a particle from site i hops to site i + 1 with rate pu{ni), and 
to site i — 1 with rate qu{ni). The parameters p and q determine the asymmetry of 
the process: The hopping is symmetric when p = q and asymmetric otherwise. For 
the function u{n) we have m(0) = 0, as particles cannot leave an empty site, otherwise 
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Figure 1. Some possible particle hopping events in the ZRP with L sites and open 
boundary conditions. 


u{n) > 0 is arbitrary. By definition the rth bond is between sites i and i + 1 for 
1 < i < L — 1. Bond 0 (L) represents the link between the lattice and a left (right) 
boundary reservoir that is not modeled explicitly. At the left boundary, particles enter 
with rate a from the left reservoir and exit with rate 7n(?7.i). Similarly, at the right 
boundary particles enter with rate 5 from the right boundary reservoir and exit with 
rate (3u{nL) (Fig. 1). 

The dynamics of the zero-range process can be conveniently represented using the 
quantum Hamiltonian formalism [9, 35]. In this approach one defines a probability 
vector \P) = En^n|n ) where | n) = |ni) 0 ... ® is the canonical basis vector 
of associated with the particle configuration n = (?7,i,?7,2,... ,ni) and the 

probability measure on the set of all such configurations. For a single site in this tensor 
product the configuration with n particles on that site is represented by the basis vector 
\n) which has component 1 at position n-|-l and zero elsewhere. By definition |P) obeys 
the normalization condition (S|P) = 1 with the summation vector 

(S| = 5^{n| (1) 


and the orthogonality condition (n|n') = (5n,n'- Within this formalism the Markovian 
time evolution of the ZRP is represented by the Master equation 

^|P(()> = -H\P[t)) (2) 

with 


H 


L-l 

ho -b hk -b flL 

k=l 

L-l 

-{Yl [P^^kht+i - dk) + flfc+i - dk+i) 

k=l 

+ a{df - 1) -b 7(ah - di) + 5(0^ - 1) + fd{dl - di)^ 


( 3 ) 




ZRP conditioned on an atypical current 


6 


where d~^ and d are infinite-dimensional particle creation and annihilation matrices 


a = 


( 0 
1 
0 
0 

V ••• 


0 

0 

1 

0 


0 0 
0 0 
0 0 
1 0 


\ 


a = 


( 0 m(1) 0 0 

0 0 m(2) 0 

000 m(3) 

0 0 0 0 

V. 




(4) 


/ 


and d = ^ diagonal matrix with the (r, s)th element given by u{r)dr^s- 

The snbscript k in and dk indicates that the respective matrix acts non-trivially on 
site k of the lattice and as nnit operator on all other sites. We also introduce the local 
particle number operator fik with the diagonal particle number matrix h = 
and note the useful identity 

y^d^y-^ = y^^d^ (5) 


for any non-zero complex number y. The global particle number operator N = 
commutes with the bulk part of H, expressing particle number conservation of the bulk 
jump processes. 

According to (2) the probability vector at time t is given by 


P(t)}=exp{-Ht)\P(t)}. 


( 6 ) 


Using (s|a+ = (s| and {s\d~ = (s|d for the one-site summation vector (s| = 
one verifies that the global summation vector (1) is a left eigenvector of H with 
zero eigenvalue. This expresses conservation of probability through the relation 0 = 
d/{dt){ S I P(t) ) = —{ S \H\ P(t)). A stationary distribution, denoted | P* ), is a right 
eigenvector of H with eigenvalue zero. 

In [29] it is shown that the stationary distribution of the ZRP is given by a product 
measure 


\p*) = \p*' 


\Pt 


(7) 


where the marginal distribution \P^) is the single-site probability vector with 
components 


P* ■=PToh[nk=n] = ^YluiP) \ (8) 

^ i=l 

The empty product is defined to be equal to 1 and Zk is the local partition function 

oo n 


n=0 i=l 

The fugacities Zk in the steady state are given by 


Zk = 


[{a + 6){p - q) - a[d +-f6] + 


L-l 


lij) - q - fd) + tdij) - q + l) j 


L-l 


( 10 ) 
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The mean density at site k is related to the fugacity by (see (8)~(9)) 



( 11 ) 


The absence of detailed balance is reflected in a steady-state current 



( 12 ) 


Notice that the existence of a steady state is determined by the radius of convergence 
of Zk. If the form of the transition rates results in some Zk outside this radius of 


convergence, then a stationary distribution does not exist. One expects relaxation to 
the local marginal distribution wherever the radius of convergence is finite, and a growing 
condensate on the other sites [29]. 

2.2. Grandcanonical conditioning 

By ergodicity, the stationary current j* is the long-time mean of the time-integrated 
current Jk(t) across a bond k,k + 1, i.e. of the number of positive particle jumps 
from site k to k + 1 minus the number of negative particle jumps from site A: -|- 1 to 
k up to time t. Because of bulk particle number conservation the stationary current 
does not depend on k. It is of great interest to consider not only the mean j*, but 
also the fluctuations of the time-averaged current jkit) = Jkit)/t, in particular its large 
deviations. The distribution gk{j, t) := Prob [jk(t) = j ] has the large deviation property 
gk{j,t) (X [33] _ Por particle systems with finite state space the large deviation 

function does not depend on the site k and satisfies the Gallavotti-Cohen symmetry 
[32, 33] and this is widely believed to be generic. However, it has been pointed out 
that neither statement remains valid for the ZRP with open boundaries if one considers 
sufficiently large atypical currents. In order to get deeper insight into this breakdown we 
study here the ZRP conditioned to produce a strong atypical current. However, instead 
on enforcing a particular (integer) value of the time-integrated current, we consider a 
“grandcanonical” conditioning defined in terms of a Legrende transform with conjugate 
parameter s. This corresponds to an ensemble of ZRP-histories where the current is 
allowed to fluctuate around some atypical mean. We refer to this ensemble as the 
s-ensemble. 

Grandcanonically conditioned Markov processes may be studied in the spirit of 
Doob’s h-transform [18, 19], as was done for the ASEP conditioned on very large 
currents in [14, 23, 24], see also [20, 21] for recent discussions. This microscopic 
approach also provides a construction to compute effective interactions for which the 
atypical behaviour becomes typical, see [14, 23] for applications to the ASEP. Besides 
the conceptual insight gained in this way into extreme behaviour, this is potentially 
interesting from a practical perspective: By adjusting interactions between particles, 
one obtains a means to make desirable, but normally rare, events frequent. 
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Following [36, 37] this construction is done by first defining a weighted generator 
H{s) where the operators that correspond to a jump to the right (left) are multiplied 
by a factor e^ Then one considers the lowest left eigenvector of H{s), i.e., the 

eigenvector (0 | to the lowest eigenvalue of H{s), which we denote by eo(s). We also 
introduce the diagonal matrix D{s) which has the components -Dn(s) of the lowest 
eigenvector on its diagonal. Since all off-diagonal elements of H{s) are non-positive we 
appeal to Perron-Frobenius theory and argue that all components of (0 | can be chosen 
to be real and non-zero. Thus D is invertible and we have 

(0| = (S|£)(s), {Q\D-\s) = {S\. (13) 

This allows us to introduce the grandcanonical Doob’s h-transform 

G(5) = D{s)H{s)D-\s) - eo(5). (14) 

Notice that by construction G{s) has non-positive off-diagonal elements and, moreover, 
by (13), (S \G{s) = 0. Hence G{s) is the generator of some Markov process, which we 
shall refer to as the effective generator, with jump rates w)n',n(s) = tCn',nT)n'('S)-Dn^('S). 
We denote by T’n(s) the invariant measure of the effective process G{s). The steady 
state of this Markov process provides the steady state of the grandcanonically biased s 
model. It is easy to prove that P^{s) = Zln(s)Fn(s) where Fn(s) are the components 
of the lowest right eigenvector of the weighted generator H{s). In order to avoid heavy 
notation we suppress in the following the dependence on s if there is no danger of 
confusion. 


2.3. Local conditioning in the ZRP 


Under the condition that the mean integrated current across some fixed bond (fc, k + 1) 
fluctuates around a certain value parameterized by s we obtain, following Ref. [37] , the 
weighted generator 

fc-l L 

H^^\s) = Y,hi + h{s)+ Y. ( 15 ) 

1=0 l=k+l 


Here 


ho{s) = - 

a(e*aj^ — 1)-|- 7 (e * 0 ]^ — di) 



hk{s) = - 

p{Gd^ dfc+i - dk) + g(e - 4+i) 

, 1 < k < L - 1 

hris) = - 

h(e — 1) + — di) 


(16) 


We refer to this setting as local conditioning. Notice that hfc(O) = hj.. We denote the 
lowest left eigenvector of by (W | = (S Specifically for A; = 0 we drop the 
subscript 0 i.e., we write {yo| = (^| = (S|Tl 

Define the partial number operator W = Yi=i From (5) one concludes 

. This implies for the left eigenvector 

{Yk\H^’^\s) = {Yk\eo (17) 
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where the lowest eigenvalue is independent of k and 

{Yk\ = {Y\e^^>^ = {S\Dk (18) 

with Dk = -De®^*. This yields effective dynamics 

= DkH^'^\s)b;^ - eo = - cq = (19) 

We conclude that the effective dynamics does not depend on k and we can focus on 
conditioning on a current across bond 0 , i.e. between the left reservoir and site 1 of the 
lattice. 


3. Microscopic density profiles and condensation 

3.1. Left eigenvector 

The left eigenvector of H^^\s) was computed in [30]. For self-containedness we repeat 
here the essential steps, which are based on a product ansatz (FI = (|/i | ® ( 1 / 21 ® ® 

{VlI = { S |T1 with the diagonal matrix D = . .yf where n are the diagonal one-site 

number operators. Then 

- ( F= ( F I I a{y,e^ - 1) + 7 ( 6 "^ - yi)dryb + 


L-l 


+ 


X] \p{yk+i - ybdkVk ^ + <i{yk - yk+i)dk+iykli 

k=l 

d{yL - 1) + - ybdiyl^ 

= ( F I I a{yie" - 1) + [7(e"" - yi) + p{y 2 - 2/i)] diyb+ 


( 20 ) 


L-l 


F Ip(!«!+i “ vO + iivt-i - Vk)] ikVt ‘ + 


( 21 ) 


k=2 


diyjL - 1) + [qidjL-i - yi) + /?(l - yr)] diyb 

One sees that {F | is a left eigenvector if the following equations are satisfied: 

0 = piyk+i - yk) + qigjk-i - yk) (22) 

0 = 7 (e“'* - yi) + p{y 2 - yi) (23) 

0 = g(2/L-i-|/ l)+/5(1-2/l). (24) 

We define the hopping asymmetry a = p/q. The ansatz 

yk = A + (25) 

yields yk-i — yk = B{a — and therefore solves (22). 

The left boundary equation (23) yields 


A — Ba^+pBa^{a — l )/7 = 0 


(26) 


which gives 


A = e 


(p - g + 7) 


Ba^ 


7 


(27) 
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On the other hand, the right boundary equation (24) yields 
B{p — q)a + /?(! — A — Ba) = 0 
Plugging this into A leads to 

— l)a~^ 


B = 


- g - /3) + /3(p - q + 
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(28) 


(29) 


7(e ‘ - l)(p -q- P) 


and 

^ ^ 7e~"(P - Q - P) + - Q + 7)q^~^ ^ ^ ^ >v- 23^7 

7(P - Q- + /^(P - Q + 'y)a^~^ 7 (p - q- (d) + l3{p - q + 7 ) 0 ^“^ ‘ 

For the lowest eigenvalue we get 

eo = - {a{yie" - 1) + d{yL - 1)) 

= -pe"{yi - 1/2) - -^qiVL-i - Vl) 

7 P 

= {p — q){—a~^e^ — ^a~^)Ba^~^^ 


(31) 


7 

from which one hnds 

eo = (p- g)(e"" - 1) 

3.2. Right eigenvector 




a(da^ ^e® — ^5 


pip - q - /d) + /d{p - q + 'y)a^-^' 


(32) 


Following [30] we make a product ansatz also for the right eigenvector: | X) = 
|(e^i) 8 ) 1^2) ® ® \xl) where \xk) is the unnormalized vector with components 


{Ql)n = xlWu{i) 


-1 


i=l 


This yields 


/ dk 

dk+i 1 


^/c+1 / 


A') 


-ft„(s)| X) = (ae‘ - 7Pi) 7"' “ ^ I IA) 

-ftilA') ={pxL-S)\l-i^^\\X) 


(33) 

(34) 

(35) 

(36) 


One sees that | X ) is a right eigenvector if the following equations are satished: 

b = pxk- qxk+i (37) 

b = ae® — ^Xi (38) 

b = (dxi — 5 (39) 

with some constant b. The ansatz 

Xk = C + Da^ (40) 
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yields pxk — qxk+i = C{p — q) and therefore solves (22) with b = C{p — q). Using the 
bonndary recursions one then obtains 


C = 


and 


D = a 


l3{p - q + 'y)a^ ^ + 7 (p - q - /3) 
ae^{p - q -/3) + 6{p - q +'y) 


/3{p - q + 7)a^“i + 7(p - q - /3)' 
Notice that in terms of the local fugacities Xk one has 

Co = —a + 7e“'^a:i — <5 + (3xl = ^(1 — e“®). 
3.3. Effective dynamics 

The transformation (5) yields for the bulk hopping terms 


(41) 


(42) 


(43) 


- DhtD-' = 




Vk 


Vk+I 


dk+1 


= p- 


Vk+i 


Vk 

+ (p - g) ( 1 


dk -dk)+q 


Vk 

'Vk+i 






fc+i 


(44) 


A 

Pk+i 


4+1 - (p - g) (1 — —) 

Vk+iJ 


dk 


and for the boundaries 

— bhQD~^ = a{yie^(x[ — 1) + 'y{ybe~^df — di) 

= ayie^dt - 1) + - di) 

+ a{yie^ - 1) + (p - g)pr^(pi - A)di 

and 

- bfiLE)-^ = diyidf - 1) + fdiybdl - di) 

= dybaf - 1) + I3yf\dl - db 

+ d(pL - 1) + (p - q)yl\yL - A)dL 
Therefore the effective dynamics is given by 
((■•"’(s) = DGl"’!)-' - Co 


(45) 


(46) 


L-1 

E 


k=l 


yk+1 


p- 


dk - 4 + g 


yk 




^/c+1 


(47) 


yk^ ^ yk+i 

ayiebdt - 1) + 7Pr^e"^(a7 - di) 

Sybdt - 1) + Pyb{di - di) 

It is intriguing that this is a driven ZRP with a spatially varying driving field Ek{s) = 
logo + 21ogi/fc+i(s)/pfc(s). This space-dependence will be present even in the case of 
non-interacting particles with u{n) = wn. Therefore, conditioning on an atypical local 
current can be realized by an effective process with a space-dependent driving field. 
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The stationary distribution of the effective dynamics is a product state 


I P-(s )) = |P,*(s)) ® \p;{s)) 0 ... 0 |P£(s)) 
where |-P^(s)) is the probability vector with components 

Tl ^ 

■ 1 
2 = 1 

Here 


( 48 ) 


(49) 


Zk = XkVk (50) 

is the local fugacity given by (25) together with (29), (30) and (40) together with 
(41) (42). The normalization Zk is the local analogue of the grand-canonical partition 
function, given by (9) as before. The density at site k is related to the fugacity as in the 
unconditioned ZRP by Eq. (11). For example, for non-interacting particles u{n) = n, 
and therefore the density is simply pk = Zk- It is remarkable that one obtains a non¬ 
trivial density prohle even for non-interacting particles. 

The stationary current is 

j*{s) = pXkVk+i-qVkXk+i = ayie" -'yxie~^ = f^xi-Syi = {p - q){AC - 


Plugging in the expressions yields 
j*[s) = {p-q) 


a(3e^a 


s^L—1 


j6e 


(3{p - q + 7 )a^ ^ ^{p - q - Id)' 


(52) 


3.4- Examples 

We now explore in more depth the results of the previous subsections for some specihc 
choices of the model parameters. 


3.4-1- Barrier-free reservoirs To somewhat reduce the parameter space of the model, 
it is natural to consider 

/3=p (53) 

7 = g, (54) 

i.e., the hopping rates out of the hrst and last site are the same as the bulk hopping 
rates. We further simplify the notation by reexpressing the parameters a and 6 in terms 
of reservoir fugacities: 

a = zi'yp/q = Zip (55) 

6 = Zr Id q/p = Zrq (56) 

where zi is the fugacity of the left reservoir and Zr is that of the right reservoir. The 
interpretation is that in the original (unconditioned) process, particles jump with the 
same rate between reservoir and boundary sites as inside the bulk of the chain. We 
therefore refer to this choice of rates as “barrier-free reservoirs”. 
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A few further notations will prove useful in what follows. First, we express the 
hopping asymmetry a = p/q via 
. L + 1 


r = 


logo. 


(67) 


It is further convenient to characterize the boundary reservoirs by chemical potentials 
Zi = e^y and to define parameters A/r, Jl and zq as follows 

A/J, = Hi - Hr = \og{zi/zr), (58) 

H ={Hi + Rr)/‘2, (59) 

Zo = e^ = (60) 

Finally, we define a relative position along the lattice 

k 


rk = 


L + 1 


(61) 


3.4-2. Partially asymmetric ZRP with barrier-free reservoirs With the choice (53)-(56) 
the parameters A, B,C,D take the simpler form 

— 1 

(62) 


A = 


^L+i _ e 


C = 


This leads to 
yk = 1-(1- 


— 1 

zie^a^~^^ — Zr 
< 7,^+1 — 1 


,L+1—fc 


B = 


D = 


— 1 
Zr — Zie^ 


a 


L+l 


sinh[z/rfc] + e ® sinh[i/(l — r^)] 




(63) 


(64) 


1 sinh [9] 

The effective hopping rates can now be read off Eq. (47). They corresponds to a ZRP 
in which the bulk bias to jump to the left and right is site dependent: a particle hops 
from site k to k -\-l with rate u{nk)pk and from site k to k — 1 with rate u{nk)qk where 


Pk=P 


yk+i rzrz sinh[z/rfc+i] + sinh[z/(l - r^+i)] 


Vk 


= ^ 


and 


Vk-x ^ 

Qk = q— = Am 


es+13 sinh[i/rfc] + sinh[z/(l — r*,)] 
es+i'sinh[z/rfc_i] + sinh[i/(l — 'rfc_i)] 


yk ' “ ■ e®+^ sinh[i/rfc] + sinh[i/(l — r^)] 

These rates correspond to a space-dependent local field 


Ek = log = 2 log 


e^+i^ sinh[z/rfc_|_i] -|- sinh[z/(l — r^+i)] 


(65) 


( 66 ) 


(67) 


qk+i “ e^+^sinhlilrt]+sinh[;^(l- r,;)] 

It is interesting to note that the effective bulk hopping rates, and thus the effective field, 
are independent of the reservoir fugacities. 

At the left boundary (site 1), according to (47), particles are effectively injected 
with rate zi Pq and removed with rate qi, where po, qi are given in (65) and (66). Similarly 
at the right boundary (site L) the effective injection and removal rates are Zr qi+i and 
Pl, respectively. Thus, the effective boundary dynamics remains of the barrier-free form 
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Figure 2. The fugacity profile Zk (left) and the effective field Ek (right) for an 
asymmetric ZRP with p = 1, q = 1/2, zi = 1, Zr = 0.4, and L = 40. Several values of 
s are plotted including the typical profile of the unconditioned process (s = 0). The 
fugacity profile and the effective field are both flat, except for boundary layers of size 
0(1) sites. 


with the same reservoir fugacities, and only the effective space-dependent hopping rates 
depend on the cnrrent-conditioning. 

We proceed to examine the typical fngacity profile dnring the atypical cnrrent event. 
To this end, we first calcnlate the right eigenvector 

— Zr sinh[z/(l — r^)] -|- sinh[z/rfc] 


Xk = 


zie^a^^^ + {zr - 


zie )a 


sinh [h] 


Thns, the fngacity prohle Zk = XkVk is 

e^ sinh^[h(l — r*,)] -f- e“^ sinh^[z/rfc] -|- 2Qsmh[i)r]^ sinh[h(l — r^)] 


Zk = Zo 

with 


sinh^ r 


Q = cosh ( s + ^ -b r 


-.( 68 ) 


(69) 


(70) 


The fngacity profile (69) along with the effective driving field are plotted in fignre 2 for 
strongly asymmetric dynamics, and in fignre 3 for weakly asymmetric dynamics. 

The cnrrent is 

. sinh(T^) sinh(z> -b s -b ^) 


j*{s) = {p-q) 


zie^a^~^^ — ZrQ * 


= 2^/^- 


-. (71) 


— 1 sinh(z/) 

This cnrrent conld be generated by a ZRP with the same constant bnlk hopping rates 
p, q as in the original (nnconditioned) dynamics bnt with effective reservoir fngacities 
z^^ = zie^^ and zf^ = ZrQ~^ ■ However, as pointed ont above, in a hnite system this 
is not how the conditioned process is typically realized. There is a space-dependent 
driving field and the effective fngacities Zk are not those of an effective left reservoir 
with bonndary fngacity zie^. 























ZRP conditioned on an atypical current 


15 



n 

• S = 1 1 
■ S = 0 

• 

♦ s = -2 

‘ • 

_ T s = -4 

__— 


Pf f f ff'fTTTT 


0.2 0.4 0.6 0.8 1 

k/(L+1) 


Figure 3. The fugacity profile Zk (left) and the effective field Ek (right) for a weakly- 
asymmetric ZRP with p = 1/2 -I- 2/L, q = 1/2, zi = 1, = 0.4, and L = 40. Several 

values of s are plotted including the typical profile of the unconditioned process (s = 0). 
The solid lines correspond to the macroscopic result, Eqs. (109) and (123). 


3.4-3. Symmetric ZRP For p = q both A and B of Eqs. (29) and (30) diverge, but yu 
is well-defined. One gets 

7(e“" - 1 ) 


1 /fc — 1 + 


vi.!^ +1) + - 1 ) 


[p + (3{L - k)] 


and similarly 


Xk = 


ae^[p + - k)] + 5[p + ^{k - 1 )] 


+ l) + - 1 ) 


For the current this yields 


]*{s) =p- 


a/3e* — yde" 


(72) 


(73) 


(74) 


vdi+1)+ha - 1 ) 

In the case of barrier-free reservoirs of /3 = 7 = p, a = zip, and 5 = ZrP, there is a 
simplihcation for which one hnds 

t/fc = e“^(l - Tfc)-h Tfc and Xk = zie^{l - rk) + ZrCk- (75) 

The fugacity prohle is therefore quadratic, 

\2 I 


A/i 2 2 , / 

Zk = zo^e ^ [1 - rk) +e ^ r^ + 2cosh{^s + — jrk{l - rk) 


zie^ - Zre ® 


and the current reduces to 

fa)=p . 

The effective held (67) in this case is 

rt+i + (1 - rt+i) 


(76) 

(77) 


El = 2 log 


e^/^ Tfc -|- e“®/2 (1 — rk) 


(78) 
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Figure 4. The fugacity profile Zk (left) and the effective field Ek (right) for a 
symmetric ZRP with p = q = 1/2, zi = 1, Zr = 0.2, and L = 40. Several values 
of s are plotted including the typical profile of the unconditioned process (s = 0) 
where the fugacity profile is linear. When s = —Ay « —1.61, the fugacity profile is 
linear again (and coincides with that of s = 0), but the corresponding effective field is 
non-zero and space-dependent. Solid lines correspond to the macroscopic result, Eq. 
(124). 


Specifically, for zi = Zr =: Zq the unconditioned ZRP is in equilibrium. Nevertheless, 
the conditioned process has non-trivial properties. There is a current j*{s) = 
2pzosmh.{s)/{L + 1) and the fugacity profile 


Zk = Zo 


1 + 4rt (1 - rt)smh^ (^ 


(79) 


depends on space in a non-linear (quadratic) fashion. The conditioned process has a 
space-dependent drift Ek = sinh(s/ 2 )[e®/^rfc-t-e“®/^ ( 1 —which decays algebraically 
with distance from the boundaries. 

The fugacity profile and effective held for a symmetric ZRP are presented in hgure 

4. 


3.4-4- Totally asymmetric ZRP Consider g = 7 = 5 = 0 where particles can jump 
only to the right. One has A = 1 and 5 = 0 (which is always true if 7 = 0). On the 
other hand, C = ae'^/p and Da’^ = 0 forl<fc<L — 1 and Da^ = ae^{p — l3)/{p/3). 
Therefore 

Xk = Zk = /p {1 < k < L — 1) 

XL = zl = aeV/9- 
For the current one has 

r{s) = ae\ ( 82 ) 

As above one may parametrize the boundary rates in terms of boundary fugacities 
a = zip and (3 = Zrp. Then Zk = zlg^ for 1 < < L — 1 and zl = zie^/zr- For the 


(80) 

(81) 
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totally asymmetric case, the effective dynamics is that of the original process, but with 
an effective reservoir fugacity zf^ = zie^. 

3.5. Spatial condensation patterns 

At this point we recall that the construction presented above is valid only for 
normalizable left and right eigenvectors, i.e., if {Y \ Pq) < oo for any given initial 
distribution |Po), (“S'| X) < oo and Zk < oo. In particular, this is the case for 
non-interacting particles where Zk is the stationary local density. On the other hand, 
depending on the choice of model parameters, the normalization condition can be 
violated. In this case, a stationary distribution does not exist and one expects the 
emergence of condensates [30, 31, 34]. It is intriguing that because of the space- 
dependence of the local fugacities the existence of such condensates will generically 
depend on the lattice site k and condensates may form in the bulk of the system. 

To illustrate how conditioning on atypical currents may lead to condensation we 
focus on the example of symmetric hopping rates, as described in section 3.4.3, and 
further restrict the case of baths with equal fugacities, zi = z^ = Zq. Assume that these 
bath fugacities are subcritical, i.e., zq < Zc, where Zc is the radius of convergence of the 
sum (9). The unconditioned process, with s = 0, has a flat fugacity profile, while any 
s ^ 0 leads to a quadratic profile (79) whose maximum is attained at site k = {L + l)/2 
and equals ^(l-i-i)/2 = 2^0 cosh^(s/2). Thus, a steady state distribution exists only if 
|s| < arccosh( 2 ^c/^o )5 and for larger |s| a condensate forms in the middle of the 
system. Similar condensation transitions occur when conditioning systems with unequal 
reservoir fugacities, or with asymmetric hopping. In all these cases, the condensates form 
in the bulk and their exact location is dictated by the maximum of the fugacity profile. 

We expect the supercritical phase to be composed of condensates within the bulk 
which continually accumulate particles, coexisting with a stationary “fluid” with a finite 
density. The exact pattern of these condensates is beyond the scope of the present paper. 

3.6. Thermodynamic limit 

Consider the limit L —)■ 00 with asymmetry a > 1 fixed and define the macroscopic space 
variable r = r^. For r 7 ^ 0,1 one obtains y{r) = 1 and z{r) = ae^/{p — q + 'j) =: z* 
independent of r. The stationary current of the effective process becomes j* = {p — q)z*. 
Therefore the effective process in the thermodynamic limit is a ZRP with constant bulk 
hopping rates p,q as the original process, but effective left injection rate a®® = ae®. 
However, there are in general boundary regions of width ^ = 1/ logo sites with spatially 
decaying local driving field and associated non-trivial exponentially decaying boundary 
layers of the fugacity. 

In the symmetric case the localization length ^ diverges. One has y{r) = r-|-e“^(l — 
r) which leads to an algebraically decaying local drift [see (78)] 

^ 4sinh(g/2) 

L[e^P r -|- (1 — r)] 


( 83 ) 
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(the L~^ factor indicates weakly asymmetric hopping, as discussed below). Therefore, 
for any s 7 ^ 0 the effective ZRP has space-dependent hopping rates. The fugacity takes 
the form 


z{r) = --K|^.^--^e V-F(l-e - r^. (84) 

1 1J 1 J 

For general s this has a quadratic term which does not exist in the ZRP with space- 
independent hopping rates. Even when the fugacity prohle is linear this is generated by 
space-dependent hopping rates. 

One may also consider the case of weakly asymmetric hopping rates, where the 
asymmetry scales as p — q = v/L. Such rates correspond to a driving field which varies 
substantially only on a macroscopic length scale. Taking the L —)■ cxd limit of Eq. (67), 
and noting that in this limit 9 ^ i>/2q, yields 

2 i/(e®’''‘" cosh[z/r] — cosh[i/(l — r)]) 

L(e^+'^ sinh[i/r] -|- sinh[z/(l — r)]) 

(here the result is presented for q = 1/2 to simplify notation). Other properties of 
weakly asymmetric systems are considered in more detail in Section 4 below. 


3.7. Time-reversed effective dynamics 


The adjoint (time-reversed) dynamics of a process with generator H is generally defined 
by H* = where the diagonal matrix P* of stationary probabilities is given 

by P* I S ) = I P*) for a stationary probability vector | P*) and transposed summation 
vector I S) = ((S |)^. The system satisfies detailed balance when H = H*. This is 
equivalent to time-reversibility and implies the absence of macroscopic currents. In 
general, however, the absence of macroscopic currents does not imply detailed balance. 
In this section we show that in our system the absence of currents does in fact imply 
detailed balance. 

Here, since P* = PT, where T has the components of the right eigenvector of H{s), 
one has 


GPs) = t{s)H^{s)t-\s)-eo{s). 


Observe that 


Therefore 


r(a+) r ^ = xPdf^, r(aj r ^ = Xkd+. 


rhfcT ^ =p[ - 4 ) + g ( - 4+1 


Xk 


^^fc+1 
Xk+1 


fXk+1 ^ , 


Xk 

+ {pXk - qXk+i) 


\ Xk 
Xk 


— g \^k^k+l 4+1 

/ Xk+1 
I ^fc+1 dk 


y^fc+l Xk 


= g- 


Xk-\-l 

Xk 


Xk 


dk^k+i dkj+P (®fc®fc+i 4+1 




( 86 ) 

(87) 


( 88 ) 
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dk+i dk 

^k+l ^k 

and for the boundaries 

— — 1) + — di) 

= ax^^e^{d'^ — di) + 7a:ie“®(a^ “ 1) + {ax^^e^ — 7)^1 + 7a:ie“® — a (89) 
= ^xie~^{di — 1) + axi^e^{di — di) + C 

and 




- f 


5{x^^d^ - 1) + (d{xLdt - dr) 

/dxiid'j^ - 1 ) + 6xl^{d1 - di) + (dxi - 5 + [dx^^ - (3)dL 

i-^^V 

XL ) 


l3xL{dl - 1 ) + Sx^^{d^ -(II) +c 


(90) 


Therefore, by the telescopic property of the sum and (43) one gets the effective adjoint 
dynamics 




L-l 



k=l 


Xk+1 

Xk 


[ 70:10 ^(aj^ 
I3xlid'l - 


(dkdk+i dkj +P (®fc®fc+i 

— 1) + Q;o:j“^e^(aj“ — di) 

1) + 5xl^{dl - di) 



(91) 


This is a driven ZRP with a spatially varying driving held = — log a + 2 \ogXk+i/Xk- 
By comparing transition rates one hnds detailed balance (j*'°^*(s) = if and only 

if the stationary current (52) vanishes. Hence, for the conditioned model, the absence 
of a macroscopic current also implies detailed balance. 


4. Atypical current events via macroscopic fluctuation theory 

It is our aim in this section to relate the exact microscopic results obtained above to the 
optimal proflles for a large deviation of the current as obtained from the macroscopic 
fluctuation theory (MFT) [5]. Since the effective dynamics does not depend on site 
k of the local conditioning we shall consider here global large deviations of the global 
time-integrated current in the whole lattice. 

The MFT can only deal with symmetric or weakly asymmetric hopping where we 
set 

p — q = v/L (92) 

and V is kept fixed in the thermodynamic limit of L —)■ cxd. In what follows we employ 
the notation defined in (58)-(60). To conform with conventional notations in studies of 
the MFT, in this section x (rather than r) denotes the macroscopic location along the 
system. 
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In the limit of L —)■ cx), and rescaling space as a; = fc/L (0 < a; < 1) and time as l/h^, 
the evolution of the ZRP may be described using a fluctuating hydrodynamic equation. 
In this description, the density p(a;, t) evolves according to the continuity equation 

d]{x,t) 


p{x,t) = 


dx 


where j{x,t) is a fluctuating current. The current is given by 

dp 


j(x,t) = -D(p)4^ +"£'((>) + 


(93) 


(94) 


where ^ is a standard white noise, and D{p), a{p) are drift and diffusion coefficients. 
The validity of the fluctuating hydrodynamics description is based on an assumption of 
local equilibrium [10], i.e., around each site k = Lx there is “box” of size 1 <C ^ T 
sites which, during times which are much larger than but much shorter than LZ (in 
the original time units, rather than the rescaled ones), is approximately at equilibrium 
with density p{x). This assumption explains why only equilibrium and linear response 
coefficients {D and a) enter the macroscopic evolution equation. 

For the weakly asymmetric ZRP dehned above with q = 1/2 and p as in (92) it is 
possible to show [7, 38] that 


« ■ 


and 


(t{p) = z{p), 

where z{p) is found by inverting the relation (11) with (9). 


(96) 


(96) 


4 . 2 . Current large deviations 


Assume that during a long period of time T, an atypical mean current j was observed 
to flow through the system. What is the probability for such an event to happen, and 
what is the most likely density prohle conditioned on such an event? These questions 
can be answered in the limit of L —)■ cx) for a general driven diffusive system which 
is described by Eqs. (93)-(94) by using the macroscopic fluctuation theory [5] and the 
additivity principle [1, 7, 10, 38, 39]. 

The probability to see an atypical mean current j during a time window T has a 
large deviations form 

P{j)^ (97) 


Assuming that the most probable way to generate this large deviations event is by a 
static density prohle p*{x), the large deviations function F{j) is given by 


m 


mm 
pO . 


IJ + 0(pP£ --xa(p)]^ 


2a(p) 


dx. 


(98) 






ZRP conditioned on an atypical current 


21 


The minimization is over all profiles which satisfy the boundary conditions dictated by 
the boundary reservoirs 

P(0) = Pi = p{zi) 

p{l) = Pr = p{zr), (99) 

where p{z) is given in Eqs. (11) and (9). The typical profile p*{x) is then precisely the 
profile for which this minimum is achieved. 

In our case, it will prove useful to describe the optimal prohle using the fugacity 
profile z{x) rather than the density. The two are related by their usual equilibrium 
relation, which for the ZRP is given in Eqs. (11) and (9). Using the fluctuation- 
dissipation relation [39] 


Idlogz 

one can transform the current large deviations function to 


F{j) = m.m.F[z{x),j] = min 

z{x) z(x) , 


2 a(z) 


dx. 


( 100 ) 


( 101 ) 


The boundary conditions now read 2;(0) = zi and z{l) = Zr- The optimal fugacity profile 
z*{x) is the one for which the minimum is achieved, and the optimal density profile is 
p*{x) = p{z*{x)). 

A general solution of the ensuing optimization problem can be obtained as follows 
[38, 39]. The functional F can be rewritten as 

F[z{x),j]= / W{z) + V{z)( — ] dx + 


'0 


where 


U(^) = 


a(z 


\dxJ 


and W{z) = 


j - ua{z') 




2 z' 


dz' — ju, (102) 


r 


8z'^ ’ ^ ^ 2 ( t ( 2 ;] 

Performing the functional derivative of F yields 


+ 


r^a{z) 


_6F _dW 
Sz dz 


dV /dz\‘^ T w N d‘^z 


dz \dxJ 

Multiplying this equation by dz/dx one obtains 
d r. /dz\‘^ 


dx - 


W{z)-V{z)(^) 


\dxJ 


= 0 , 


which upon integrating yields the differential equation 



W{z) + K 


1 + 


iya{z) 

J 


+ Ka{z). 


(103) 


(104) 


(105) 


(106) 


V{z) cr(z) ^ 

Here K and K are integration constants which are related by iP = j‘^K/2 + jv. They 
are determined by the boundary condition ^(1) = Zr. Using (106) one may simplify a 
bit expression (102) evaluated at the optimal prohle: 


F{]) = 2 / W{z*)dx + 


2U 


J \ .'2 ) 

-dz' + 


ra[z'] , , j^K 


(107) 
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Substituting the drift and diffusion coefficients (95) and (96) in (106) yields the ordinary 
differential equation 


dz , ^ , 


i+—y+Kz. 




(108) 


As in the microscopic calculation, the fugacity profile is independent of the rates u{n), 
while the density profile depends on the form of u{n) through (11) and (9). 

The solution of (108) with the boundary conditions 2 ;( 0 ) = zi and z{l) = z^ is 

e^ sinh^[z/(l — x)] + e“^ sinh^[i/x] + 2 Qsinh[ra] sinh[z/(l — a:)] 


Z*{x) = Zq 
with 


sinh^ V 


Q = \ 1 + 


j sinh z/ 2 
ZqU 


(109) 


( 110 ) 


When integrating Eq. (108), the boundary condition 2 ;( 1 ) = Zr was used to find that 


K = ^ 

i 


1 + 


ZqU 


Q coth V 


cosh( A/i/2) j 


( 111 ) 


j sinh v\ sinh v 

An interesting limit is that of a symmetric ZRP, in which z/ —)■ 0. Taking this limit 
in Eq. (109) yields the simpler quadratic expression 

z1^q{x) = Zo[e^{l - xf + e~^x^ + 2 a;(l - x)^Jl + {]/ZqY]. ( 112 ) 

This expression simplifies further in the equilibrium case where zi = Zr, i.e., A/i = 0: 




=0^) = *o[l + 2x(l - n(\/l + (iAci)^ - !)]■ 


(113) 


Note that according to (109) the fugacity profile (and hence the density profile) is 
symmetric around a: = 1/2 whenever zi = Zr, or equivalently A/i = 0, even when the 
bulk dynamics is asymmetric (i.e., when z/ 7 ^ 0 ). 


4-4- Optimal profile in the s-ensemble 

As discussed above, rather than considering a system conditioned on a specific value of 
the current, one may move to an ensemble in which the current is allowed to fluctuate but 
is “tilted” towards this atypical value. Formally, this is done by considering the cumulant 
generating function of the current,(i.e., its log-Laplace transform) (^(s), defined as 

qTG{s) ^ _ j 

where (97) was used. The last integral is dominated, in the large T limit, by its saddle 
point, and therefore G{s) is the Legendre transform 

G'(s) = max[sj - F{j)]. (115) 

J 

In this grandcanonically conditioned “s-ensemble”, the optimal profile for a given 
value of s is exactly the optimal profile for the value j(s) at which the maximum in Eq. 
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(115) is attained (this is nothing bnt the eqnivalence of the two ensembles). This valne 
can be fonnd by inverting the relation 

s{j) = F\j). (116) 

We shall now carry ont this compntation for the ZRP. Note that the cnrrent j enters 
the optimal prohle (109) only throngh Q given by (110), and therefore obtaining Q{s) 
is goal of the cnrrent calcnlation. 

The first step is to find F{j), by snbstitnting the optimal prohle in (107). Carrying 
ont the integrals, one hnds 
rl 


z*{x)dx = 


zo 


and 


u sinh u 

dz 


Jo z*{x) 3 

For the ZRP, one also obtains 


coth r — 1) + cosh( A/i/2) (cosh z/ 


= -log[(5+ \/Q‘^ - 1]. 


sinh V 


J - Tza[z 
2 z' 


■dz' = 


jdz' v{zi - Zr) 

2 z' 2 


- Zr) j 


Zl 


Combining Eqs. (103), (107), (111) and (117)-(119) yields 


(117) 

(118) 

(119) 

( 120 ) 


m = 


ZQl' 


cosh log( - 1 + Q) .(121) 


sinh u 

Then, nsing (116), we hnd 

ZqI' . , / A/^ 

sinh s H- 


j(s) = 


sinh r 




+ r 


( 122 ) 


(123) 


Finally, snbstitnting in (110) yields the simple expression 

Q = cosh ^ + iJj , 

which together with (109) yields the optimal fngacity prohle for a given valne of s. As 
before, it is interesting to consider the symmetric limit of u = 0. In this case, 

zI^q(x) = zo[e ^(1 — x)^ + e ^ x‘^ + 2 coshes + ^J^x{l — x)]. (124) 


4-5. Comparison with the microscopic calculation 

As discnssed in Section 3.6, in the weakly asymmetric limit of (92) with q = 1/2, and 
when L —)• oo, one has z/ —)• z/ [where z/ is the asymmetry parameter dehned in (57)]. 
Therefore, in this limit the microscopic fngacity prohle (69)-(70) is precisely the same 
as the macroscopic expressions (109) and (123). Similarly, the macroscopic expression 
(122) for the cnrrent in the s-ensemble coincides in this limit with exact expression (71) 
obtained for barrier-free reservoir conpling with weak asymmetry. Indeed, hgnres 3 and 
4 present fngacity prohles for symmetric and weakly-asymmetric hopping and several 
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values of s, and demonstrate the agreement between the microscopic and macroscopic 
results. 

An interesting observation is that the microscopic expression (69) for the optimal 
profile and macroscopic one (109) have the same functional form even when L is small 
and the asymmetry is not weak, provided one uses the definition x = k/{L + 1). In 
particular, by substituting u = cL in the macroscopic prohle (109) (with a constant c) 
it is possible to obtain the correct profile for strongly asymmetric hopping (i.e., with 
p ^ q independent of L). However, using the macroscopic theory beyond its limit of 
validity comes with a price: The precise value of the prefactor c which corresponds to 
given values of p and q cannot be determined within the macroscopic approach, and 
only the microscopic one yields equation (57). 
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